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We present a new independent evaluation of the hadronic and QCD contributions to the QED running coupling 
a(Mz) and to the muonium hyperfine splitting v. We obtain: Aahad = 2770(17) x 10~ 5 and A^had = 232.5(2.5) Hz. 
Combined with the QED and Electroweak Standard Model contributions, they lead to: a~ 1 (Mz) = 128.926(25) 
and to the Fermi energy splitting i/p — 4 459 031 783(229) Hz , where for the latter, we have used, in addition, the 
precise measurement of the muonium hyperfine splitting ^ CX p. We use vf in order to predict the ratios of masses 
m^/m e = 206.768 276(11) and of the magnetic moments [i^j [i% = 4.841 970 47(25) x 10 -3 , which are in excellent 
agreement with the ones quoted by the Particle Data Group. These remarkable agreements can provide strong 
constraints on some contributions beyond the Standard Model. 

1 Introduction 

A natural and important extension of our recent work Q (hereafter referred as SN) on the hadronic and QCD 
contributions to the anomalous magnetic moment of the muon and tau leptons ai , is the evaluation of the hadronic 
and QCD contributions to the QED running coupling a(Mz) and to the muonium hyperfine splitting v. The 
analysis is important as it will give a complete set of estimates of three independent observables within the same 
inputs used in SN. These self-contained results should pass tests from a comparison with high precision experiments 
and different existing predictions in the literature. Using a dispersion relation, it is remarkable to notice that the 
different lowest order hadronic contributions for these three processes can be expressed in a closed form as a 
convolution of the e + e~ — > hadrons cross-section au{t) with a QED kernel function K{t) which depends on each 



observable: 
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where a 1 (0) = 137.036 and Mz = 91.3 GeV. It behaves for large t like a constant. 
— For the muonium hyperfine splitting Ai^h a <j, the kernel function is (see e.g fl|): 



where: 




and we take (for the moment) for a closed comparison with Q| [], the value of the Fermi energy splitting: 

VF = 445 903 192 0.(511)(34) Hz . (8) 

It behaves for large t as: 
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The different asymptotic behaviours of these kernel functions will influence on the relative weights of different 
regions contributions in the evaluation of the above integrals. 

2 Input and Numerical Strategy 

The different data input and QCD parametrizations of the cross-section a nit) have been discussed in details in 
SN jij and quoted in the last column of Table 1, corresponding to the estimate in different regions. Table 1 is 
analogous to Table 2 of SN. We shall only sketched briefly the numerical strategy here: 

— Our result from the 1 = 1 isovector channel below 3 GeV 2 is the mean value of the one using r-decay and e + e~ 
data. In both cases, we have used standard trapezoidal rules and/or least square fits of the data in order to avoid 
theoretical model dependence parametrization of the pion form factor. In the region 0.6 — 0.8 GeV 2 around the 
uj-p mixing, we use in both cases e + e~ data in order to take properly the SU(2)p mixing. The SU(2) breaking in 
the remaining regions are taken into account by making the average of the two results from r-decay and e + e~ and 
by adding into the errors the distance between this mean central value with the one from each data. 

— For the 1 = isoscalar channel below 3 GeV 2 , we use the contributions of the resonances w(782) and 0(1020) 
using narrow width approximation (NWA). We add to these contributions, the sum of the exclusive channels from 
0.66 to 1.93 GeV 2 . Above 1.93 GeV 2 , we include the contributions of the w(1.42), w(1.65) and 0(1.68) using a 
Breit-Wigner form of the resonances. 

— For the heavy quarkonia, we include the contributions of known J/ip (IS to 4.415) and T (IS to 11.02) families 
and use a NWA. We have added the effect of the tt bound state using the leptonic width of (12.5 ± 1.5) keV given 
in§]. 

— Away from thresholds, we use perturbative QCD plus negligible quark and gluon condensate contributions, 
which is expected to give a good parametrization of the cross-section. These different expressions are given in SN. 
However, as the relative role of the QCD continuum is important in the estimate of Aahad, we have added, to 
the usual Schwinger interpolating factor at order a s for describing the heavy quark spectral function, the known 
a 2 mg/t corrections given in SN. However, in the region we are working, these corrections are tiny. 

— On the Z-mass, the integral for Aahad has a pole, such that this contribution has been separated in this case 
from the QCD continuum. Its value comes from the Cauchy principal value of the integral. 

3 The QED Running Coupling a(M z ) 

The result given in Table 1 corresponds to the lowest order vacuum polarization. Radiative corrections to this 
result can be taken by adding the effects of the radiative modes 7r°7, r]j,ir + ir~j, .... We estimate such effects to 
be: 

Aa had = (6.4 ±2.7) x lO" 5 (10) 

by taking the largest range spanned by the two estimates in (5) and Q. Adding this (relatively small) number to 
the result in Table 1, gives the total hadronic contributions: 

A« had = 2769.8(16.7) x 10~ 5 . (11) 
Using the QED contribution to three-loops || : 

Aa QED = 3149.7687 x 10" 5 , (12) 
and the Renormalization Group Evolution of the QED coupling: 

a' 1 (M z ) ^«- 1 (0)[l- Aa QED - Aa Uad ] , (13) 

the next section, we shall extract this value from the analysis. 
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Table 1: Lowest order determinations of Aahad x 10 and Ai/h a d [Hz] using combined e + e and inclusive r decay 
data (2nd and 4th columns) and averaged e + e _ data (3rd and 5th columns). 

Region in GeV 2 Aai ia d x 10 5 A^had [Hz] Data input 



Light Isovector 

Ami ~> 0-8 
0.8 -► 2.1 
2.1 -> 3. 
Total Light 1=1 

Average 
Light Isoscalar 
Below 1.93 



0.66 -> 1.93 
From 1.93 to 3 
w(1.42), w(1.65) 
0(1.68) 

Total Light L=0 
Heavy Isoscalar 

J/ip(lS^> 4.415) 
T(15-> 11.020) 
T(349) 

TbiaZ Heavy 1=0 
QCD continuum 

3. -> 4.57 2 
4.57 2 -> 11.27 2 
11.27 2 -> (M z -3 
(M z + 3 GeV) 2 -> 
Z-pole 
4M t 2 -> oo 
ToiaZ QCD Cont 



r+e + e 

314.5 ±2.3 

77.2 ±3.4 

62.3 ±9.2 



302.7 ±7.1 
82.0 ±5.4 
53.6 ±4.9 



T+e + e 

152.9 ± 1.8 
12.1 ±0.5 
7.8 ± 1.2 



148.4 ±3.1 
16.9 ± 1.9 
6.7 ±0.6 



454.0 ± 10.6 438.2 ±10.2 112.8 ±2.2 172.1 ±3.7 



GeV) 2 
4M t 2 



446.1 ± 10.4 ± 7.9 



31.5 ±1.1 
52.3 ± 1.2 

11.6 ±3.0 

9.4 ± 1.4 

14.6 ±4.6 

119.0 ±5.9 

116.3 ±6.2 

12.7 ±0.5 
-(0.1 ±0.0) 
128.9 ±6.2 

330.1 ± 1.0 
503.0 ±1.0 
2025.7 ±2.0 

-(794.6 ±0.6) 

29.2 ± .5 
-(24.0 ±0.1) 
2069.4 ±5.2 



172.5 ± 3.0 ± 0.3 



12.7 ±0.4 
13.7±0.3 
2.7±0.7 

1.2 ±0.2 
1.7 ± .5 

32.1 ±1.0 

4.0 ±0.2 
0.1 ±0.0 

W 

4.1 ± .2 

17.5 ±.1 
5.0 ± .1 

1.3 ±0.0 

W 
W 
W 

23.8 ± 1.4 



[3, 7, 

\l \\ 
ill 



NWA |] 
NWA || 
^exclusive [|10f 




NWA 
NWA 
NWA 



[91 
|9| 

111 



(u, d, s) 
(it, d, s, c) 
(it, d, s, c, 6) 

principal value 

(«, d, s, c, 6, i) 



Final value 



2763.4 ± 16.5 



232.5 ±3.2 



one obtains the final estimate: 

aT 1 (M z ) = 128.926(25) , (14) 

which we show in Fig 1 for a comparison with recent existing determinations. One can notice an improved accuracy 
of the different recent determinations j|, [| 0, [ll| 0, which are in fair agreement with each others. Also a detailed 
comparison of each region of energy with the most recent work of || shows the same features (agreement and slight 
difference) like in the case of a M in SN, due to the slight difference in the parametrization of the data and spectral 
function. However, the final results are comparable. Finally, one can remark that due to the high-energy constant 
behaviour of the QED kernel function in this case, the low-energy region is no longer dominating. For a a, the 
contribution of the p-meson below 1 GeV is 68% of the total contribution, while the sum of the QCD continuum 
is only 7.4% (see e.g. SN). Here the situation is almost reversed: the contribution of the p-meson below 1 GeV is 
only 2%, while the sum of the QCD-continuum is 73.6%. For this reason, improvement due to the new Novosibirsk 
e + e~ data jL5) in the low-energy region will not be very significant. At present, new BES data |L6| in the J/tp 
region are also available, which can be alternatively used. Below the J/tp resonances, the BES data are in excellent 
agreement with the QCD parametrization to order a^. used here for 3 flavours, justifying the accuracy of your 
input. Above the J/ip resonances, the parametrization used here (sum of narrow resonances ±QCD continuum 
away from thresholds) can also be compared with these data. On can notice that, in the resonance regions, the 
BES data are more accurate than previous ones, which may indicate that our quoted errors in Table 1 for the J/tp 

2 Previous works are quoted in |l 4[- 
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YT01 

SN 01 
(this work) 




2700 



2800 



2900 



128.85 128.90 



128.95 129.00 



Figure 1: Recent determinations of Aahad and a (Mz)- The dashed vertical line is the mean central value. 
References to the authors are in || 0, fLj|. 



family contributions are overestimated. In addition, the threshold of the QCD continuum which we have taken 
above the J/ip resonances, matches quite well with the one indicated by the BES data. We expect that with this 
new improved estimate of a(Mz), strong constraint on the Higgs mass can be derived. 



4 The Muonium Hyperfine Splitting 



Our final result from Table 1: 



Ai/ ha d = (232.5 ± 3.2) Hz 



(15) 



is shown in Table 2 in comparison with other determinations, where there is an excellent agreement with the most 
recent determination M. Here, due to the (ln£)/£ behaviour of the kernel function, the contribution of the low- 



Table 2: Recent determinations of Ai/h 



ad 



Authors 


Ai/ had [Hz] 


FKM 99 Q 


240 ± 7 


CEK 01 ||| 


233 ±3 


SN 01 (This work) 


232.5 ±3.2 



energy region is dominant. However, the p-meson region contribution below 1 GeV is 47% compared with the 68% 
in the case of a M , while the QCD continuum is about 10% compared to 7.4% for a M . The accuracy of our result is 
mainly due to the use of the r-decay data, explaining the similar accuracy of our final result with the one in Q using 
new Novosibirsk data. The agreement with (J] canbe understood from the agreement of the averaged correlated 
e + e~ and r-decay data compiled in || with the new Novosibirsk data used in Q. We differ with DH98 |ll| in 
the treatment of the QCD contribution 3 . For light quarks, QCD is applied in the region where non-perturbative 
contributions are inessential. For heavy quarks, perturbative QCD is applied far from heavy quark thresholds, 
where it can be unambiguously used. Adding to this result, the QED contribution up to fourth order, the lowest 
order electroweak contribution n , and an estimate of the higher order weak and hadronic contributions [fL8| : 



A^qed = 4 270 819(220) Hz 

Avyjc^tl n) = — m„m... I I 



{1.0) = 



-65 Hz , 



3 For more details, sec 
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Ai/bad(M ~ 7(2) Hz , (16) 
one obtains the Standard Model (SM) prediction: 

fSM = vf + Ai/q ED + At/ wcak + Aj/ had + Az/ had (/i.o) . (17) 

If one uses the relation: 

VF = PF ( I 7— 7 , (18) 

with: 

16 

p F = — {Za) 2 Z 2 cR 00 , (19) 

one would obtain: 

v>SM = 4 463 302 913(511)(34)(220) Hz , (20) 

where the two first errors are due to the one of the Fermi splitting energy. The first largest one being induced by 
the one of the ratio of the magnetic moments. The third error is due to the 4th order QED contribution where, 
one should notice that, unlike the case of a^, the dominant errors come from the QED calculation which should 
then be improved. Unfortunately, these previous errors are still too large and obscure the effects of the electroweak 
and hadronic contributions. One the opposite, the data are very precise |19[]: 



^cxp 



= 4 463 302 776(51) Hz . (21) 



Therefore, at present, we find, it is more informative to extract the Fermi splitting energy vp from a comparison of 
the Standard Model (SM) prediction with the experimental value of v. Noting that Vp enters as an overall factor 
in the theoretical contributions, one can rescale the previous values and predict the ratio: 

— = 1.000 957 83(5) . (22) 

Vp 

Combining this result with the previous experimental value of v, one can deduce the SM prediction: 

v F = 4 459 031 783(226) Hz , (23) 

where the error is dominated here by the QED contribution at fourth order. However, this result is a factor two 
more precise than the determination in Q given in Eq. (^), where the main error in Eq. (|^) comes from the input 
values of the magnetic moment ratios. Using this result in Eq. (E3h into the expression: 



m e \ 1_ 

mj (1 + m e /m M ) E 



v F = PF\zr-) T, , ~ 7ZI N 3 (1 + °m) , (24) 



where: 

p F = ^{Za) 2 Z 2 cR OQ , (25) 

and Z = 1 for muonium, a- 1 (0)=137.035 999 58(52) |o), cR^ =3 289 841 960 368(25) kHz |l) and a M = 
1.165 920 3(15) x 10~ 3 (^2|, one can extract a value of the ratio of the muon over the electron mass: 

— ^ = 206.768 276(11) , (26) 
m e 

to be compared with the PDG value 206.768 266(13) using the masses in MeV units, and with the one from Q: 
206.768 276(24). After inserting the previous value of m e /m^ into the alternative (equivalent) relation: 

VF=PF ^ 



one can deduce the ratio of magnetic moments 

Ms 



^- = 4.841 970 47(25) x 10~ 3 , (28) 



compared to the one obtained from the PDG values of /x^//i p and [i p j pf B Mm/Mb = 4.841 970 87(14) x 10~ 3 . 
In both applications, the results in Eqs. ( p6| ) and ( p8| ) are in excellent agreement with the PDG values. These 
remarkable agreements can give strong constraints to some contributions beyond the Standard Model and should 
be explored. 
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5 Conclusions 



We have evaluated the hadronic and QCD contributions Aahad and A^had respectively to the QED running coupling 
and to the Muonium hyperfine splitting. Our results shown in Eqs. (|ll|) and ([H]), are in excellent agreement with 
existing determinations shown in Fig. 1 and Table 2 and are quite accurate. These results have been obtained 
within the same strategy and data input as the one of the anomalous magnetic moment obtained previously in SN 
[Q. For this reason, they are self-contained outputs. One of the immediate consequences of these results is the 
prediction of a(Mz) given in Eq. ([l4|), while we have used the result for the muonium hyperfine splitting for a 
high precision measurement of the ratios of the muon over the electron mass given in Eq. (|26| ) and of magnetic 
moments given in Eq. (|2S|) . These Standard Model predictions are in excellent agreement with the ones quoted 
by PDG |9|]. These agreements can be used for providing strong constraints on some model buildings beyond the 
Standard Model. 
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